function ca=delta(x,y)
s=size(x);
for i=1:(s(2)-1)
    dx=x(i+1)-x(i);
    dy=y(i+1)-y(i);
    ca(i)=cos(atan(abs(dy/dx)));
end
ca=[ca 1] ;


